Resonant Superfluidity in an Optical Lattice 



I TitvinidzeS M Snoek^ and W Hofstetter^ 

^ Institut fiir Theoretische Physik, Johann Wolfgang Goethe-Universitat, 60438 
Frankfurt am Main, Germany 

^ Institute for Theoretical Physics, Valckenierstraat 65, 1018 XE Amsterdam, The 
Netherlands 

E-mail: ireikliSitp . uni-f rankf urt . de 
FAGS numbers: 37.10.Jk, 67.85.Fq, 67.85.-d 

Abstract. We study a system of ultracold fermionic Fotassium (^°K) atoms in a 
three-dimensional optical lattice in the vicinity of an s-wave Feshbach resonance. 
Close to resonance, the system is described by a multi-band Bose-Fermi Hubbard 
Hamiltonian. We derive an effective lowest-band Hamiltonian in which the effect of 
the higher bands is incorporated by a sclf-consistcnt mean-field approximation. The 
resulting model is solved by means of Generalized Dynamical Mean-Field Theory. In 
addition to the BEG/BCS crossover we find a phase transition to a fermionic Mott 
insulator at half filling, induced by the repulsive fermionic background scattering 
length. We also calculate the critical temperature of the BEC/BCS-state and find 
it to be minimal at resonance. 
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1. Introduction 

The first experimental realizations of Bose-Einstein condensates in dilute atomic gases 
of rubidium lithium [2] and sodium [3] atoms initiated a new field of condensed- 
matter research, providing an ideal laboratory for comparing theoretical models and 
experimental results with high accuracy. In particular the important consequences of 
Bose-Einstein condensation could be investigated, which up to 1995 had remained an 
elusive and inaccessible phenomenon in experiments. 

Not long after the first realization of BEC, ultracold fermionic gases were studied 
experimentally as well. The first important results of quantum degeneracy in trapped 
Fermi gases were obtained in 1999 at JILA [1] and later on by other groups O [6]. A 
break-through experiment in this field was the investigation of fermionic superfluidity 
at the crossover between the BEC state and the Bardeen-Cooper-Schrieffer (BCS) state 
[3 El El Uni [11]. This is made possible by the use of Feshbach resonances, which have 
become an indispensable experimental tool for ultracold atom experiments. Feshbach 
resonances not only allow one to tune the interatomic interaction with high precision, 
but also make it possible to increase it to a level where the critical temperature becomes 
high enough for the investigation of attraction-induced superfluidity. On the contrary, 
away from resonance the critical temperature is usually exponentially suppressed and 
experimentally inaccessible. The regime of strong, even diverging interactions, the 
so-called unitarity region, on the other hand defines a new field of research where 
standard mean-field methods break down and the physics has to be described in a non- 
perturbative way. Moreover, this system allows for the experimental investigation of 
the BEC-BCS crossover: for negative scattering lengths the system is a BCS superfiuid, 
whereas for positive scattering length fermionic atoms with opposite spin pair up to 
form a bosonic molecular bound state. More recent experimental work has focused on 
studying the effect of spin imbalance on the BCS state, i.e. the case when an unequal 
number of fermions occupies the two different spin states [121 1131 [El [IS], as well as 
mixtures of fermions with unequal masses, such as ^Li and [T7] . 

Also the effect of periodic potentials on trapped Fermi gases has been studied 
experimentally [IHl [IHl 1201 [21] • Recently, evidence for a fermionic Mott insulator 
was obtained in a system of repulsively interacting fermions in an optical lattice 
[22l|23]. Optical lattices and Feshbach resonances have been combined experimentally 
as well: the group at ETH Ziirich reported the production of molecules in three- 
dimensional cubic optical lattices using s-wave Feshbach resonances in early 2006 [19], 
but no evidence of a superfiuid state in the lattice was found until later that year, 
when superfiuid ^Li was loaded in an optical lattice at MIT where both a condensate 
and an insulating state were observed [20]. The results where interpreted in terms of 
a superfiuid-to-Mott insulator transition, for which a detailed theoretical description is 
still lacking. 

In this work we study an ultracold mixture of fermionic atoms in two different 
hyperfine states in a three-dimensional optical lattice close to a Feshbach resonance. 
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Figure 1. Schematic picture of the BEC-BCS crossover in an optical lattice. By 
tuning the interaction strength between the two fermionic spin states, one obtains a 
smooth crossover from a BEC regime of tightly bound bosonic molecules (a) to a BCS 
regime of large Cooper pairs (c). In between these two extremes, one encounters an 
intermediate crossover regime where the pair size is comparable to the interparticle 
spacing (b). For total fermionic filling one, the system can undergo a phase transition 
to the Mott insulator phase (d). 



This system has all the characteristics of the continuum BEC/BCS crossover described 
above: For magnetic field values below the resonance, fermions with different spin form 
bosonic molecules (see Fig. [T]). By varying the magnetic field the bosonic level is detuned 
relative to the fermionic one, which changes the ratio of the densities of fermions and 
molecular bosons as well as the the effective interaction between the fermions. On top 
of this BEC-BCS crossover physics, which is familiar from the system without lattice, 
new features emerge when an optical lattice is applied. The most prominent one is 
the occurrence of a fermionic Mott insulator for half-filled fermions deep in the BCS 
regime, which is stabilized by the repulsive fermionic background scattering length. As 
described below, we find a first order transition between the BEC/BCS state and the 
Mott insulator. For a total filling of two fermions per site the BCS state competes 
against a band insulating state [2H |25l [26] . 

The presence of an optical lattice allows to utilize powerful non-perturbative 
methods that are available for lattice systems. Here we apply Generalized Dynamical 
Mean-Field Theory (GDMFT) [271 I2H] to the system described above. However, 
GDMFT is a single-band approach, whereas Feshbach-resonant interactions in an optical 
lattice lead to a multi-band model |3ll [351 133 [371 EH |39] . We therefore first perform 
a mean-field decoupling of the higher bands, thereby deriving an effective single-band 
Hamiltonian which is self-consistently coupled to the higher bands. 

The paper is organized as follows: in the next section we introduce the microscopic 
model and in Sec. |3]we introduce the Generalized Dynamical Mean- Field approach we 
use to solve this model. In Sec. |4| we present the result of our numerical calculations 
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and in Sec. [5] we end up with concluding remarks. In the appendix, we describe in detail 
how the self-energy for the resonantly interacting Bose-Fermi mixture studied here can 
be calculated in the Dynamical Mean-Field framework. 



2. Microscopic Model 

Studying ultracold fermions close to a Feshbach resonance is a challenging problem. 
Due to the fact that exactly on resonance the scattering length is infinite, the standard 
fermionic Hubbard Hamiltonian cannot be defined. To solve this problem, it is necessary 
to formulate a two-channel Hamiltonian by separating out the resonance state and 
treating it explicitly [10]. The nonresonant contributions give rise to a background 
scattering length. As the Feshbach resonance occurs due to a coupling with the bosonic 
molecular state, the additional degrees of freedom introduced in the two-channel theory 
are bosons [lO] . 

An ultracold atomic gas of fermionic atoms and molecular bosons close to a 
Feshbach resonance in the presence of an optical lattice is thus well described by a Bose- 
Fermi Hubbard model [351 SI] • Iii our calculation we assume the molecular bosons to be 
in the lowest band. For the fermions, on the other hand, we have to take into account 
also the higher bands, in order to properly incorporate the two-body physics associated 
with the Feshbach resonance [3Sl ESI 112] • Since the bandwidth is much smaller than the 
band gap, we approximate the higher bands to be flat and only take into account the 
full band-structure for the lowest band. Moreover, we neglect the interaction between 
fermions in higher bands with each other and with the bosons. This is justified because 
the filling in the higher bands is very small, so that interaction effects are also small. 
The Hamiltonian thus has the following form: 

oo 

n =h5 + h, + h;, + 5^(h^ + h^,), (1) 

1=1 

^'f = -^fH 4,os-..o + UfY. KtAifi - (/^ - ^) E ' (2) 

(ij) i i 



(ij) 



^% = Ufb ^Ho + 9oJ2 (^1^,^,0^4,0 + h.c) , (4) 

i i 

i 

where c]^ i is the creation operator of a fermion with spin a for the /-th band on lattice 
site i. h\ is the creation operator of a boson at site i. h{^i = cl^iC^fj-i, nf^ = fii^^i + hii^i are 
the fermionic number operators, and = 6-6^ is the bosonic number operator, if and 
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tb are the tunneling amplitude for fermions and bosons, respectively, f//, Uh and Ufh are 
Fermi-, Bose-, and Bose-Fermi-Hubbard interactions, respectively. These interactions 
arise due to the background scattering lengths. Furthermore is the chemical potential, 
5 is the detuning of the bosonic level, and u is the frequency of the harmonic oscillator 

associated with an optical lattice well. Finally, gi = go \J Lp''^-'(0) is the Feshbach 
coupling to the /-th band of the lattice, where go is the Feshbach coupling for the 
lowest Hubbard band and L^^^'^\q) is the generalized Laguerre polynomial. Choosing 
the Feshbach couplings in this way guarantees that the two-body physics associated 
with the Feshbach resonance is incorporated exactly [35l US] . 

The parameters of the generalized Hubbard Hamiltonian ([T]) are given by: 

Vo 




exp 




(7) 



rrih/mf 



3/4 



V 27rh) 



Vo\ 

E'rJ 



3/4 



(8) 
(9) 

(10) 



Here e/^^^ 



= A/iinag(-B - 

h"^ / 2X^171 f(b) is the recoil energy, Vq is the amplitude of the optical lattice 
potential and A is the laser wavelength, a/, a^, and afb are the fermion-fermion, 
boson-boson, and fermion-boson background scattering lengths. In our calculation we 
approximate the background boson-boson and the Bose-Fermi scattering lengths by 
ttfe = 0.6a/ [43] and a/6 = 1.2a/ Furthermore, B is the magnetic field, and Bq 

and AB are the position of the Feshbach resonance and its width, respectively. A/Xmag 
is the difference in the magnetic moment between the closed and open channel of the 
Feshbach resonance. Finally, m/ and nib are the respective masses of the fermions and 
bosons. 

The Hamiltonian ([T]) can be simplified by the following rescaling: 

Shu , , 

/^ = /^-^, (12) 



5 = 5- 



2 



(13) 



such that the factor ^ disappears. 



3. Method 



3.1. Derivation oj the effective single-band Hamiltonian 



The multi-band Hamiltonian derived so far is very complicated, since it both involves 
strong correlations and many bands. Simply neglecting the higher bands would lead to 
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an incorrect description close to the Feshbach resonance, since the Feshbach parameter 
g is very large and even exceeds the band gap [35]. However, the filling of fermions in 
the higher bands is strongly suppressed by the band gap. This allows us to perform a 
mean- field decoupling in the higher bands [35]. The lowest band is left untouched in 
this procedure, since the fermionic filling can be large there. 

We thus perform the following decoupling for Z > on each site: 



njli 



(14) 



This step implies that the lowest band and the higher bands are only coupled in a mean- 
field way. They can thus be diagonalized separately, but are coupled by the mean-field 
self-consistency relations. The full Hamiltonian is now given by 



Li 



where the following terms are added to the bosonic part of the lowest band Hamiltonian: 



1=1 



(16) 



where the mean-field A has been defined as A = —Yl'iLidiA^i'^ixi)- each of the 
higher bands / > we obtain the following Hamiltonian (here we suppress the site index 



2lf\w — ^ ~9i{b) 
-gi{b^) -{mu-fi) 




(17) 



The system described by Eqs. (15), (16) and (17) needs to be solved self-consistently 

'a ' ' ' 1—1 

with respect to the mean fields (6) and A. 

To diagonalize the Hamiltonian (17), one has to perform the following Bogoliubov 
transformation: 

2lhijj — ft —gi{b) 

-gi{b^) -{2lhw-fi) 






where 



u, 



UlVi 



i2lhuj - /i)2 + gf 



1 2lhuj - 

2 ^ 2^1 
1 



2 

9i(b) 
2ui 



2lhjjj — fi 

2uJi 



This leads to the following expectation values: 



n, 



2vf + 2iu]-vf)fiui) 



2lhuj — fi 



tanh 



\2kTJ 



9i{b) 
2ui 



tanh 



2kT) 



(19) 
(20) 

(21) 
(22) 

(23) 
(24) 
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where f{uJi) is the Fermi function and T is the temperature. We use absolute values in 
the equation for (ci^cii) because of the ambiguity of the sign, arising from the fact that 
still a divergence has to be subtracted (see below). 
The total number of fermions is equal to: 



,f:(:--^ta..(^)). (25) 



= "-0 

1=1 

This is a converging sum, which can be evaluated numerically. 
From Eq. ( [l6| follows that we have to evaluate the sum 

^9,(oA> = ±WE|;'-"(2^)- (26) 

1=1 1=1 

which is divergent. This divergence always arises in the gap equation of the BCS model 
when the T-matrix is approximated by a delta-potential [35l H2] . One way to solve 
this problem is by using a pseudo-potential instead of the delta-potential [42j. Here, 
however, we follow Ref. [55] and explicitly isolate the diverging contribution from the 
sum. 

First, we notice that for large /, coi can be approximated by ui = 2lhu — Ji and 
tanh (^) — 1. Therefore 

9 / 9 OO 9 

y^tanhf^)^(y^tanhf^)+ V , 

^2ui \2kT) \^2ui \2kT) ^ 2(2lhuj - ft) 

1=1 ' \i=l ' l=N+l ^ ^' 




^''''^(2fcT) X^2(2Znt-/i) +£2(2//^-/.)) ■ ^^^^ 

Here is a large integer number (in our calculation we took = 500). 

The first two terms of Eq. (27) are finite sums, but the last term diverges. This 
sum is known from the literature ^21 H5] . To separate the diverging part, we have to 
take the following limit: 

V 9\ _^ 9lLr\^) 9lL'^"\r) 
^2{2lhio-fi) ^2{2lhuj-Ji) r^o ^ 2{2lhuj - ft) 

= - lim (dv^n-fi/2n.)/r(-p'/L - 1/2) _ VI ^ ps) 



Since the diverging part ^ jr is independent of the model parameters, we can cure the 
divergence by neglecting this term [351 132] • Doing so, we obtain 

'gl^V[-Jilhuj)lV{-Jilhuj - 1/2) 



1=1 ^ 

N 2 ^ 



2hu 



^ ^ .-y^tanhf^^) . (29) 



^2{2lhuj-fi) f-^2ui \2kTJ 



We now fix the sign by requiring A > 0, since this solution minimizes the (free) energy. 
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Summarizing, we have reduced the multi-band problem to an effective single-band 
Hamiltonian: 



Tt — ''-/ ~r "T rtf, -\- Tt J 



(30) 



where 7^°, l-ib, "Hjf, and are given by Eqs. (^2), (3), (4h and (16), respectively. 



The chemical potential /i has to be adjusted such that the total filling is equal to 
the desired value ntot: 

2lhuj 



2n^ + 



(=1 



— tanh ( — ) 
\2kT) 



ntot 



(31) 



This leads to the following self-consistency loop: we start from an initial guess of 
the superfiuid order parameter (b) and calculate A using Eq. (29). As a result we know 



all parameters in the Hamiltonian (30), and can find its eigenvalues and eigenvectors. 



and correspondingly calculate new correlation functions, including the superfiuid order 
parameter (b). With this step the self-consistency loop is closed. 

It is worth noting that the effective single-band Hamiltonian we have derived here 
is different from the effective single-band model in terms of dressed particles derived in 
other approaches [361 EH]: the bosons and fermions in our Hamiltonian correspond to 
the bare particles in the lowest band. 



3.2. Generalized Dynamical Mean-Field Theory 



To analyze the Hamiltonian (30) we use Generalize Dynamical Mean Field Theory 
(GDMFT) which is explained in detail in Ref. [271 I2H]- Here we only mention that 



within GDMFT one considers a single site which is self-consistently coupled to a 
dynamical fermionic bath corresponding to DMFT [291 130] and a static bosonic mean- 
field corresponding to bosonic Gutzwiller [211 E21 E3]- These are the leading order 
contributions in a 1/ ^-expansion of the effective action {z being the lattice coordination 
number). Hence GMDFT is exact in infinite dimensions and expected to be a good 
approximation for the cubic lattice considered here, for which z = 6. The typical 
accuracy for low-temperature expectation values is around 20 percent. 

In the specific case considered here the system is described by a generalized single 
impurity Anderson model (GSIAM) with the following Anderson Hamiltonian: 



And 



H 



-TjAnd 



njAnd < njAiia < nj 



/And 



And 



(32) 



H n (n 

2 ^ 



(2/i - 6)h'' 



H 



And 



{ztb(p + A)6"^ + h.c 
U fhh^h^ + Qq {b\c^c^ + h.c^ 

UfHnl + J2{^kal^ka + [cla^^ + h.c) } + W^fc ( 



^k'\^k^ ~1~ h.c. 



where z is the lattice coordination number and = (6) is the superfiuid order parameter. 
k labels the noninteracting orbitals of the effective bath, Vk are the corresponding 
fermionic hybridization matrix elements, describes the superfiuid properties of the 
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bath and a^^ is the creation operator of a fermion in the k^^ orbital of the bath with 
spin a. Ufa = c\^a is the fermionic number operator and h = hf^^ + fif^. 

To solve the Anderson Hamiltonian we use exact diagonalization as impurity solver 



[HI HZl HH HI] . In this algorithm the infinite number of orbitals in the Hamiltonian (32) 
is truncated and only a finite number of Ug orbitals is considered. The resulting finite- 
size problem is fundamentally different from the finite-size problem of a finite number of 
lattice sites of the original Hubbard model, and the truncation procedure can be viewed 
as using a finite number of parameters (energy scales) to describe the local dynamics 
encoded in the Weiss Green's function: 

ns _|_ £- 

^And(«^n) = 'SaXndi^^n) =tUJn + Jl + Y^ ^2 i " 2 i U/2 ' ^^^^ 

where (3 is the inverse temperature and u;„ = (2n + l)7r//3 are the Matsubara frequencies. 

To close the self-consistency loop by using the lattice Dyson equation we calculate 
the normal and superfiuid Green's functions which can be written as follows 



C* — e 

G{iUn) = G^{iujn) = I deD{e)— , ^ ' (^5) 

-oo \s ^\ + ^SC 



oo 



FiiLOn) = -T.sc{iu^n) I deD{e)- ^ , (36) 

J-oo 1^ ~ ^1 + ^SC 

where C = i^n + fi — ^{ioJn) and D{e) is the non-interacting density of the states of the 
cubic lattice. S(zci;„) and T.sci'i^n) are the normal and superfiuid self-energies, which 
as shown in Appendix A can be expressed via a set of higher order Green's functions: 

T.{iUn) = T.„{i0Jn\ 



{UfQffa{i!^n) + UfbQfba{i(^n) + (J gQ*g^„{iUn)) Gl{iUn) 



G^{iuJn)Gl{iuJn) + F{aiujn)F*{aiujn) 
{aUfQff^aa{iuJn) + (^UfbQfbaa{iuJn) + gQ*g^{iujn)) F*{aiLJn) 
^ G^{icOn)GUiuJn) + F{aiuJn)F*{aiuJn) ' ^ ' 

^ {UfQff^{iUn) + UfbQfbr{iUJn) + gQ*gl^{i(^n)) F{i(^n) 



G^{iUn)Gl{iUn) + F{iUn)F*{-iUn) 
[UfQfmii^n) + UfbQfbuii^n) + gQ*g^{iujn)) G^{iUn) 



G^{iuJn)Gl{tUJn) + F{iuJn)F* 



(38) 



Here G{iujn) = {{c^fl, cI^q))^) and F{iujn) = {{c^ q-,Cio))ui are the normal and superfiuid 
single particle Green's functions. In addition we have also defined the following 
additional interacting Green's functions: Qffa{iuJn) = {{faflfa^fl))^, Qffaa{i(^n) = 

{{Iblfi))^ and Qg^.itiOn) = {{fjb^h)).- Here 

((^' B))^ = -y Y.(^\A\m){m\B\n)- ± — , (39) 

Z ^ Em- En- lUJn 
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and 



(40) 



is the partition function. 

The relation between the Weiss field and the Green's function is given by the local 
Dyson equation: 

g-^lUJn) = S(^C^n) + G-\lUn) , (41) 

where 

G{iUJn) F{iuJn) 

F{lUn) -G*{lUJn) 

is the matrix of interacting Green's functions, 



G{iujn) 



S(iw„) = 
is the self-energy matrix and 



(42) 




(43) 



(44) 



is the matrix of Weiss Green's functions. 

To determine new parameters for the Anderson Hamiltonian, we fit the Weiss 



functions calculated from (33) and (34) to the ones calculated from the eigenstates 



of the Anderson Hamiltonian via the local Dyson equation (41). We use a steepest 
decent method with the following norm: 



X 



2(-/V'max H" 



^ ' max 



n=0 



(45) 



where A^max is the number of Matsubara frequencies taken into account. 

The minimization procedure works as follows: we start from an initial guess of 
the GSIAM parameters (e^o-, Via and VF/), and then knowing the local Green's functions 



calculated from Eq. (39) and the self-energies calculated from (37) and (38) we calculate 



the lattice Green's function according to Eqs. (35) and (36). Subsequently, using 
the Dyson equation (41) we can calculate the Weiss Green's functions Q~lx{iuJn) and 
J-'~l^{icOn) ■ The next step is to fit this result by the parameterization in Eqs. (33) and 



(34) and thus to find a new set of parameters for the GSIAM. These new parameters 



serve as input for the next iteration. This procedure is repeated until convergence is 
reached. 
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3. 3. Calculation of the critical temperature 

The combination of the mean-field approximation in the higher bands and GDMFT 
explained so far, however, leads to a problem. Both approximations involve the 
superfluid order parameter The mean- field approximation for the higher bands 

implies that the local correlator (&lcj^^;Cjj^^/) is approximated by (Cj^^/Cj^^^^)- The 
GDMFT scheme, on the other hand, involves the approximation to replace the non- 
local correlator (})\hj) by {bl){bj). This means that (b) both measures the local phase 
coherence between bosons and fermions and the non-local bosonic long range order. 
However, these are two very different quantities which generally cannot be described 
by a single mean-field order parameter. At zero temperature, this problem is not too 
severe, because in this case one expects both long-range order and on-site boson-fermion 
coherence, so that (6) is large for both reasons. At finite temperature, however, this 
becomes a real problem, because the bosonic long range order is expected to vanish at 
temperatures of the order of the bosonic hopping tb- The local boson-fermion coherence, 
on the other hand, persists for much higher temperatures, since the coupling g is orders 
of magnitudes larger than tf,. Indeed, we find that in the approximation outlined above 
the full GDMFT calculations are in good agreement with a single site approximation. 
In the single-site approach the impurity site couples neither to the fermionic nor to the 
bosonic bath and long range order cannot be inferred. The critical temperature obtained 
from this calculation can be identified with the pair breaking temperature Tpair, which 
is much higher than the relevant temperature in experiments. 

The scheme explained in the previous subsection can therefore not be used to infer 
the critical temperature for superfluid long range order. In order to do so, we have to 
modify the approximation and remove the ambiguous nature of the order parameter 
(b). This is made possible by the observation that the term Ab^ in the Hamiltonian 
merely renormalizes the self-energy of the bosons: in the BEG regime the bosons are in 
a coherent state and this term is equivalent to a shift of the bosonic chemical potential. 
This is also clear from the treatment in [33], where terms from higher bands enter the 
bosonic self-energy. To make this more explicit, we write 

A = - ^gM,c,,) = ±{b) ( — 

1=1 ^ 

+ y^J^ ,-y^tanhf^)U(6)A'. (46) 

We can therefore replace the term 

- (^AP + /i.e.) = - {A!{b)b^ + /i.e.) (47) 
in the Hamiltonian, by 

- A'¥b , (48) 

such that terms from the higher bands only renormalize the chemical potential. We 
remark here, that this might look like an additional approximation. However, one has 
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to keep in mind that the term connecting A to the bosonic creation operator originated 
from the mean-field approximation in the higher bands. By treating the contribution 
of higher bands within the bosonic self-energy we are therefore restoring part of the 
mean-field approximation made in the previous step. Indeed, by using second order 
perturbation theory in the couplings to higher bands (which is justified if the band 
energy exceeds the Feshbach coupling), we can also obtain this correction directly as 
part of the bosonic self-energy, without invoking a mean-field decoupling. 

This improved approximation for treating higher bands gives for T = similar 
results as before; in particular the position of the transition to the Mott insulator is in 
good approximation the same. The superfiuid order parameter is smaller, as expected. 
However, for nonzero temperatures this improved approximation scheme allows for a 
calculation of the critical temperature for superfiuid long range order, which was not 
possible in the previous approximation. 
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Figure 2. Bosonic and fermionic filling (a) and superfiuid order parameters (b) as a 
function of magnetic field B for T = 0. The dotted line corresponds to the Feshbach 
resonance, while the dashed line indicates the phase transition from the superfiuid 
phase into the Mott insulator phase. The magnetic field is measured in units of Gauss. 



4. Results 

We study a mixture of potassium atoms (^°K) and Feshbach molecules in a three- 
dimensional optical lattice. The on-site harmonic oscillator frequency is chosen to be 
u = 271 X 58275Hz, which corresponds to a lattice with wavelength A = 806nm and 
Rabi frequency of Qr = 27r x 1.43GHz. The Feshbach resonance considered here is 
at -B = 202. IG and the width of the resonance is 7.8 G [50j. The difference between 
the magnetic moments of the closed and open channels of the Feshbach resonance is 
A/i = 16/9/2^, where hb is the Bohr magneton. The total filling per lattice site in our 
calculation is ntot = 1- 
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4.1. Zero temperature 

First we consider the case of zero temperature. Our calculations for the ground state 
are summarized in Fig. [2j Deep in the BEG regime only bosonic molecules are present. 
When the magnetic field is increased, close to resonance the number of fermions is 
increasing and the number of bosons decreasing. Above the resonance we mainly have 
fermions and the number of bosons is small. Both fermions and bosons are superfluid. 
We remark again that here we describe the physics in terms of bare bosons and fermions: 
in terms of dressed particles as in Ref. [35], these are still molecular bosons and the 
BEC/BCS crossover takes place when the bosonic self-energy crosses twice the Fermi 
energy [35]. However, in the case of half- filled fermions, this crossover is intercepted 
by a first order phase transition to a fermionic Mott insulator state, which happens 
at a critical value of the magnetic field of -B = 249G. Calculations which include only 
the lowest band of the Bose-Fermi-Hubbard model (as well as with one and two exited 
bands), yield this transition into the Mott insulator phase already close to the Feshbach 
resonance at i? ^ 205G (results not shown). This implies that to capture the superfiuid 
region 205T ^ B < 249T higher bands which renormalize the bosonic self-energy are 
crucial. 

Another point worth noting is the first-order nature of the transition to the Mott- 
insulating state. In contrast, if one integrates out the bosonic degree of freedom and 
describes the Feshbach resonance in terms of an effective, attractive interaction between 
the fermions within a single channel model for the lowest band, one finds a different 
scenario. In this case the induced attractive interaction dominates, until it is cancelled 
by the repulsive background interaction. This means that within the single channel 
approximation one finds a regime with a normal Fermi-liquid phase in between the 
BEC/BCS phase and the Mott insulator, which is absent in our phase diagram. This 
directly indicates that the effect of higher bands is crucial to capture the first-order 
transition between the superfluid and insulator phase. 

4-2. Nonzero temperature 

Having clarifled the ground state phase diagram, we now consider flnite temperature. 
In particular we investigate the critical temperature for the transition to the normal 
state. Deep in the BEC regime, the critical temperature is constant (T^ ~ 0.21tf) and 
completely determined by the properties of the bosons: the bosonic hopping parameter 
tb, the interbosonic background scattering length and the bosonic density. Only very 
close to resonance the critical temperature suddenly drops (see Fig. |3]). This coincides 
with the magnetic field value for which fermions enter the system. On the BCS side 
of the resonance, the critical temperature depends on the magnetic field and increases 
with B (see Fig. [3]). This implies that at resonance the critical temperature is minimal. 
This is in sharp contrast to the situation where no lattice is present, in which case the 
critical temperature is maximal close to resonance. 

This surprising fact can be understood from the behavior of the critical temperature 



Resonant Superfluidity in an Optical Lattice 



14 




0.05 0.1 0.15 0.2 0.25 150 160 170 180 190 200 210 220 230 240 



T B 

(a) (b) 

Figure 3. Finite temperature results: in subfigure (a) we plot the fermionic superfluid 
order parameter as a function of temperature T for different magnetic fields above the 
resonance, while in subfigure (b) we show the phase diagram. The blue solid line 
separates the superfluid phase from the normal phase. Here temperature is measured 
in units of the fermionic hopping tf. 



in the single-band attractive Hubbard model [STl [52] . In this model, the critical 
temperature is low both for very large and very small attraction and has a maximum 
in between. The reason for the low critical temperature at small attraction is the 
conventional exponential suppression of in the BCS regime. For strong attraction 
the critical temperature decreases again because the fermions start forming bound pairs 
with a greatly enhanced effective mass. Identifying the resonance position with the 
case of very large attraction, this explains the low critical temperature at this point. 
When moving away from resonance, the effective attraction induced by the Feshbach 
resonance becomes weaker and hence the critical temperature increases again. Far 
away from resonance one would therefore expect to find a maximum of the critical 
temperature, after which it decreases again because the BCS regime of weak attraction is 
entered. However, due to the transition into the Mott insulator phase for the unit filling 
considered here we cannot see the maximum of the critical temperature. Estimating the 
induced attractive interaction at the transition point to the Mott insulator by assuming 
it to be equal to the repulsive background interaction, this indeed gives a value for the 
induced attractive interaction which is larger (in absolute value) than the position of 
the maximum in the single-band attractive Hubbard model [5T1 [52] . 

Our calculation also shows that on both sides of the resonance the ratio [c^C]) / (b^) 
as a function of the temperature for fixed value of the magnetic field is constant. This 
means that the on-site Bose-Fermi coherence is not affected by the temperature for the 
low temperatures considered here. 
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5. Summary 

We have studied ultracold fermionic atoms in a three-dimensional optical lattice 
close to a Feshbach resonance. We derived an effective description in terms of a Bose- 
Fermi Hubbard model, in which the molecular degree of freedom is explicitly present. 
Our calculations show in agreement with Ref. [35] that the effect of higher bands is 
crucial for a correct description of the Feshbach physics. We therefore take into account 
the fermionic occupation of higher bands. 

To solve the strongly interacting multi-band problem we decouple the higher bands 
from the lowest one via a mean-field decoupling and reduce the Hamiltonian to an 
effective single-band Bose-Fermi Hubbard model which is self-consistently coupled to 
the higher bands. To solve this resulting model we use GDMFT. 

The low temperature physics close to the Feshbach resonance is very rich. Upon 
changing of the magnetic field, the ratio of fermionic and bosonic densities is changing. 
Below the resonance the system is mainly occupied by molecular bosons forming a 
condensate. Close to resonance the number of bosons decreases while the number 
of fermions increases. The fermions are in the superfiuid phase. This resembles the 
BEC/BCS crossover close to a Feshbach resonance without an optical lattice. In 
addition, for the unit total filling considered here we found a transition into the fermionic 
Mott insulating phase when the magnetic field is increased even further. The Mott 
insulator phase is stabilized by the repulsive fermionic background scattering, which 
at large magnetic fields overcomes the attractive interaction induced by the Feshbach 
resonance. The phase transition into the Mott insulator is found to be of first order. 
We found that higher bands are crucial for a quantitatively correct prediction of the 
transition point. 

We also calculated the critical temperature of the BEC/BCS superfiuid phase 
across the resonance. Below resonance the critical temperature is independent of the 
magnetic field, until it sharply drops close to resonance. Above the resonance the critical 
temperature is increasing again, leading to the remarkable result of a minimal critical 
temperature at resonance. 
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Appendix A. Derivation of the self-energy for the Bose-Fermi mixture 

In this appendix we evaluate the self-energy via correlation functions. For this purpose 
we use the equation of motion, which in general has the following form: 



B)h = { 



A,B 



) 



(A.i; 



Here Un are the Matsubara frequencies and ((A, B))^^ is the general form of the Green's 
function, with the usual notation [A, B]± = AB ± BA] the plus sign applies when both 
operators are fermionic, otherwise the minus sign is used. 

For the resonantly interacting Bose-Fermi mixture the generalized single impurity 
Anderson Hamiltonian has the following form: 

^And ^ _ J^fif.ni + Uffi^nl + Uftnffi' + g (/f/ffe + h.c}j + n^'"' 

fccr fccr k 

where fl and are the fermionic creation operators on the "impurity site" and 
in the band respectively. ¥ is the bosonic creation operator on the impurity site, 
n-^ = n| -|- n| = flfai ""-^ — ii^'^'^ is the bosonic part of the Hamiltonian. 

To calculate the self-energy we first evaluate the following commutator relations 







k 


(A.3) 


^And^ /t" 


_ = -l^fJl + Uf fifth + UfJlPl 


' - (^gfab'^ + VkaCl^ , 
k 


(A.4) 




= -SkaCka - Vkafa " <^Wkcl- , 




(A.5) 


-V/And 


= ^kaci^ + Vkafl + (rWkC^- , 




(A.6) 



where a 



-a. 



Now we use the equation of motion (A.I) for the case when A = fa and B = f^. 



In combination with the commutation relation (A.3) we get: 

{zOOn + /i/.) {{fa, fl)). - Uf{{fJlf-a. fl)). 

- UMfihJl))^ + ag{{flhJl)). -Y.^ka{{tkaJl)). 



(A.7) 



To calculate {{cka, w)uj we again use the e quation of motion Eq. (A.I), but in this case 

ation: 



with A = Cka and B = f^. With Eq. (A.5) we obtain the following re^ 

(Za;„ - Ska) {{Cka, fl)h - Vka{{L Pa)) ^ " ^Wk{{cl, fl)) 



0. 



Finally to calculate {{clg., fl))uj we use Eq. (A.I) with A = c|.- and B = /J, which 
results in 



(za;„ + Sk,) {{cl, fl))^ + Vk-a{{fl fl)). - aWk{{c,^, fl)). 



0. 



(A.9) 
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From Eqs. (A. 8) and (A. 9) we derive 



aVkaWk 

{iUn - Ska){iUJn + Ska) " 



{{fai fcr))iJ 



(A. 10) 



Now we combine Eq. (A. 10) with Eq. (A. 7) and obtain: 

{iun + fifa - A„{iun)) G„{iujn) - Asc{(^i^^n)F* {-aiUn) 

- UfQffa{i^^n) - UfbQfba{iuJn) " crgQ*^^{iUn) = 1 . (A.ir 

Note that {{fa^fl))w = Gu{iujn) is the normal Green's function and {{f^,f^))uj = F'i^] 
the superfluid Green's function. We also define 

iUn + £ka 



12) 

{eka-iU}n){£ka + iU}n) + W^' 



A 



A*sc{-iUr, 



E 



2 ' 



(A.13) 



as the normal and the superfluid hybridization functions respectively, and the following 



correlation functions: Qjfaii'^n 



{{fafafa, fa))LO, Qffa 



{{fafo-fay fa))i- 



QfbaiiUJn) = {{f^b^bJl))^, QfbaS^n) = ( /a) QgaiiUJn) = {{f^b^fl))^ and 
Qgaa{iUJn) = {{Lb\fa))cv 

To obtain the self-energy we need to derive one more equation. For this purpose, 
we again use the equation of motion Eq. (A.l) and take A = and B = f^. Based on 
Eq. ( |A.4[ ) we get: 

{tur. - lij^) ((/t, fl))^ + Uf{CflflL fl)). + UMifM fl))u - crg{{Uh\ fl))^ 

+ Y.'^kMlJl)). = ^- (A.14) 



We now replace ((c|.^, fl))^ using Eq. (A. 10) and obtain: 

- a {iun - fifa + A* F*{(TiUn) + crAsc{-<7i^^n)G^{iUn) -UfQ 

- UfbQ}i,aA^uJn) - crgQg^{iUn) = . 



(A.15) 



We proceed to write our results in matrix form. For a = 1 we use Eq. (A.ll) and 
the complex conjugate of Eq. (A.15), while for a = — 1 we take Eq. (A.15) and the 
complex conjugate of Eq. (A.ll): 

{iun + /i/t - A^{iuj„)) G^{iuJn) - AsciiuJn)F*{-iuJn) 

-UfQff^{iujn) - UfbQfb-riiuJn) - gQ*gi^{iuJn) = 1 , 
- {iun - /i/4, + A^(-iw„)) GKiUn) - Asc{iuJn)F{iuJn) 

- UfQ*fj-^{iuJn) - UfbQ*fblii^n) + gQgui^i^n) = 1 , 

{iUn + /U/t - A^{iUn)) F{iun) + Asc{i(^n)G\{iUn) 

- UfQff^^l{iujn) - U fbQ fbui^^^n) - gQ*g^{iujn) = , 
{iun - fJ^fi + A^,{iun)) F*{iun) - Asc{iu}n)G^{iUn) 

-UfQ)j^^^{iuJn) -UfbQ*fbl't('^^n) + gQg^ii^n) = . 



Resonant Superfluidity in an Optical Lattice 18 

The last four equations can be rewritten in matrix form in the following way: 
1 \ ^ / iUJn + fJ^f^ - A^{iUn) -Asc{iu}n) \ f G^{iuJn) F{iujn) \ ,^ 

UfQff^iiuJn) + UfbQfhr{i(^n) + gQ*gi^{iujn) UfQff^^iiiUn) + U fbQ fbui^ujn) + gQ*g<^{iu}n) 



Now we compare Eq. (A. 16) with the Dyson equation, which has the matrix form: 

g'\lUn) - ±{tCOn) = G-\lCOn) , (A.17) 

where 

is the matrix interacting Green's function, 

^^^^^ = ( + - At(^c.„) -Asc{^u^n) \ ' ^^^^g) 

is the matrix Weiss Green's function and S(i^) is the matrix self-energy. From this 
comparison it follows directly that 



UfQff^iiUn) + UfbQfb^{iu}n) + gQ*g<^^{iu}n) UfQff^^iiiUn) + UfbQfbni'^(^n) + gQ*g^{iujn) 

^fQ)f,i^i^^n) + UfbQ*fbi^{iU}n) - gQcn{i(^n) UfQ)f^{iuJn) + U fbQ)bi{i^n) " gQgui^^n) 



. -1 

G^{iujn) F{iuJn) 
From here we obtain the final result: 



(A. 



[UfQffaiii^n) + U fbQ fbaii^^n) + (rgQ*g^„{iUn)) G^ 



lUJr. 



Ga{iuJn)G%{iuJn) + F{aiuJn)F*{aiuJn) 
^ [aUfQff^„^{iUn) + (TUfbQfbaS^n) + gQlS^n)) F*{aiuJn) (A 21) 
G„{iUn)G%{iujn) + F{aiujri)F*{aiujri) 



(UfQff^iiUn) + UfbQfbt{'>'^n) + gQ*gl^{iUJn)) F{iuJn) 

^""^'"^"^ ~ G't(za;„)G*(zu;„) + FitUn)F*{-tuJn) 

{UfQff,ui'^^n) + U fbQ fbui^^^n) + gQ*g:^{iU}n)) G^{iuJn) 

G^{iUn)Gl{iUn) + F{iujn)F*{-iujn) 
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